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ABSTRACT 

A wide range of recent observations have shown that AGN-driven cavities may 
provide the energy source that balances the cooling observed in the centres of "cool 
core" galaxy clusters. One tool for better understanding the physics of these cavi- 
ties is their observed morphological evolution, which is dependent on such poorly- 
understood properties as the turbulent density field and the impact of magnetic fields. 
Here we combine numerical simulations that include subgrid turbulence and software 
that produces synthetic X-ray observations to examine the evolution of X-ray cavities 
in the absence of magnetic fields. Our results reveal an anisotropic size evolution of 
that is dramatically different from simplified, analytical predictions. These differences 
highlight some of the key issues that must be accurately quantified when studying 
AGN-driven cavities, and help to explain why the inferred pV energy in these regions 
appears to be correlated with their distance from the cluster center. Interpreting X- 
ray observations will require detailed modeling of effects including mass-entrainment, 
distortion by drag forces, and projection. Current limitations do not allow a discrim- 
ination between purely hydrodynamic and magnetically-dominated models for X-ray 
cavities. 
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1 INTRODUCTION 



Feedback from active galactic nuclei (AGN) is widely be- 
lieved to be the most promising mechanism to solve a 
number of problems in our understanding of the history 
of baryonic structure formation, includin g: (i) re^ulatin 
I supermassive black hole formation (Chura zov et al.l l200 
JSiiacki Springell l2006h : (ii) setting the upper mass cut- 
off in the galaxy ma ss function, and determining its evo - 
lution with redshift (iBenson et all l2003l : iBest et al.ll2005l l: 
(iii) providing the nongravitational "preheating" that is 
likely to have been important in determining the entropy 
profi l es of non "cool-core" gal axy clusters (jPonahue et al.l 
I2OO6I : IScannapieco 6hll2004l l: and (iv) oflFsetting the ra- 
diative losses seen in cool- core galaxy clusters (see e.g., 
iMcNamara NulsenI |2007[ ). Such clusters show strong 
peaks in their central X-ray surface brightness distributions, 
indicating that the central gas is cooling r apidly. Yet the 
deficit of star formation and <1 keV gas fe.g.. lRaflFertv et aP 
I2OO6I ) means that radiative cooling must be balanced by 
an unknown energy source. Currently, the most successful 
model for achieving this balance relies on heating from a cen- 
tral AGN; yet the details of this process are poorly under- 



stood (e.g ., "Brugsen k KaiseJ I2OO2I : [Reynolds enil I2OO2I : 
iBrighenti Mathews 2006). 

While AGN in cool-core clusters are observed 
to drive large bubbles into the intracluster medium 
jMcNamara e t al 



(ICM) 



200d: 



. .,B lanton et al.1 I2OOII : 

Fino guenov JonesI l200l[ iNulsen et al.1 1200^ ). the syn- 
chrotron radiation emitted by the relativistic electrons in 
these bubbles fades and becomes extremely diflficult to 
detect after about 10^ years. Moreover, the corresponding 
X-ray surface brightness depressions are only visible near 
the centre of the cluster where the contrast is large. AGN 
have also been observed to induce shocks and/or sonic 
motions in the IC M that are believed to dissipate their 
energy int o this gas (Fabian et al. 2003; Fabian et al. 2 0051: 
Kraft et al. 2 007; Ruszkpwski e t ah 2004; McNama r a et*^ 
20051 : iFinoguenov et all I2OO8I : ISanders FabianI I2OO8I I. 



although the impact of the resulting heating is difficult 
to quantify observationally. Thus, it is unclear how far 
AGN-driven cavities rise in the cluster, how they couple to 
the surrounding medium, and how they evolve. 

More specifically, the presence of these cavities has 
raised a number of key questions. The buoyant bubbles in- 
flated by the central AGN are unstable to the Rayleigh- 
Taylor (RT) instability, which occurs whenever a fluid is ac- 
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celerated or supported against gravity by a fluid of lower 
density. Yet, these cavities appear to be intact even af- 
ter inferred ages of sever al 10^ yrs, as is t he case for the 
outer cavities in Perseus (jNulsen et al.|[2005 ). On the other 
hand, purely hydrodynamic simulations fail to reproduce 
these observations. Instead, the RT and other instabilities 
shred the bubbles i n a re l atively short time (Briigge n et aL 
I2OO5I : iHeinz et all I2OO6I : IPizzolato Sokeil l2006h . Mag- 
netic fields have been shown t o alleviate this problem 
somewhat J^obinson et al. 200i; |j ones De Youn j l2005l : 
iRuszkowski et al.ll2007f ). but they also reduce the extent to 
which the interior of the hot bubbles couples to the sur- 
rounding medium, making it much more difficult for AGN 
heating to balance cooling. 

In a recent paper, Scannapieco & Briiggen (2008) 
showed that although pure-hydro simulations indicate that 
AGN bubbles are disrupted into pockets of underdense gas, 
more detailed modeling of turbulence indicates that this is 
a poor approximation to a cascade of structures that con- 
tinues far below current resolution limits. Using a subgrid 
turbulence model developed by Dimonte & Tipton (2006), 
they carried out a series of simulations of AGN heating in 
a cool-core cluster with the adaptive mesh refinement code, 
FLASH. These simulations showed that Ray leigh- Taylor in- 
stabilities act on subgrid scales to effectively mix the heated 
AGN bubbles with the ICM, while at the same time pre- 
serving them as coherent structures. The AGN bubbles are 
thus transformed into hot clouds of mixed material as they 
move outwards in the hydrostatic medium, much as large 
airbursts lead to a distinctive "mushroom cloud" structure 
as they rise in the hydrostatic atmosphere of Earth. This 
allows X-ray cavities to remain intact out to large distances 
from the cluster centre while still coupling to the surround- 
ing medium. 

Alternatively, it has been suggested that instead of un- 
derdense pockets of ideal gas, the cavities are produced by 
magnetically domina ted jets ( Li et al. 2006). In 3D MHD 
simulations by Nakamura et al. (2006, 2007) bubbles are in- 
ffated by a current-carrying jet that injects magnetic ffux 
into a small volume in the vicinity of the supermassive 
black hole. The jets are launched by injecting non-force- 
free poloidal and toroidal magnetic ffelds, and form large 
currents which travel along the inner jet axis into the lobes 
and return on the outer boundary of the lobes, forming a 
sheath around the jet axis. Like bubbles of hot, underdense 
gas, such jets expand subsonically into wide lobes that ap- 
pear cooler than the surrounding medium. Unlike hot, un- 
derdense gas, these magnetically dominated bubbles behave 
differently as they rise thr o ugh t he cluster atmosphere. 

Recently, "Piehl et al.1 (|200 8h investigated whether one 
can use the measured sizes of X-ray cavities observed at dif- 
ferent locations in their host clusters to discriminate between 
these two models. Compiling the sizes and radial offsets of 64 
cavities in the X-ray halos of clusters and groups, they were 
able to show a tight correlation between these two quan- 
tities, which is substantially different than one would ex- 
pect from simple analytic estimates in the pure-hydro case. 
From this comparison, they came to the preliminary con- 
clusion that the data favor the current-dominated magneto- 
hydrodynamic jet model. 

While a useful first step, such analytical prescriptions 
neglect some important physical effects that are likely to 



affect the interpretation of the data. While the loss of pres- 
sure due to bremsstrahlung radiation is negligible, mass en- 
trainment produced by hydrodynamic instabilities is likely 
to be important (^Pavlovski et al.i »2007i) . This will add to 
the growth of the bubbles produced by the expansion dur- 
ing the ascent in the stratified atmosphere. In fact, there is 
some circumst antial evidence for m ass entrainment in FR I 
radio sources. ICroston et al.l (|2008h showed that equiparti- 
tion internal pressures are typically lower than the external 
pressures acting on the radio lobes, so that additional non- 
radiating particles must be present. A correlation between 
the structure of the radio sources and the apparent pres- 
sure imbalance can be taken as observational evidence that 
entrainment may provide this missing pressure. Moreover, 
drag forces and acceleration by buoyancy distort the bub- 
bles, complicating the interpretation of bubble radii. 

Related to these issues is the question as to how to best 
measure the energetics of AGN-infiated bubbles. It is now 
standard procedure to use the size of a bubble and the am- 
bient pressure at its location to get an estimate of the pdV 
work needed to infiate the bubble, which is then used to i nfer 
the e nergy associated with the outburst (e.g., Birzan et al.l 
1 2 0041 ) . Interestingly, if one plots the inferred energy as a 
function of radius for the known X-ray cavities, one finds 
that the energy increa ses as a function of bubble distance 
from the cluster centre (Rafferty et al. 2006, see Fig.[8|). This 
effect is not only apparent in a large sample of different clus- 
ters, but even within the few individual clusters with known 
multi ple bubble s, such as Hydra A and Perseus. Recently, 
Sand ers Fab ian (2007) reported the possible detection of 
another cavity in Perseus at very large radii (~ 170 kpc), 
and inferred an outburst energy several times higher than 
that inferred for bubbles at ~ 30 kpc from the cluster cen- 
tre. Even larger differences between inner and outer bubbles 
have been reported in Abell 2204 where the outer bubbles 
have radii of 240 kpc and 160 kpc (jSanders et al.ll2008l ). Un- 
derstanding this effect will be crucial in correctly addressing 
the question whether AGN can be a general solution to the 
cooling ffow problem in clusters and galaxies. 

Here we combine adaptive-mesh refinement (AMR) 
simulation s that include a subgrid turbulence model 
(Scann apieco &: Briiggenl I2OO8I ) with software to produce 
synthetic X-ray observations to better r eproduce the de- 
tailed evolution of hydrodynamic bubbles (IHeinz Briiggenl 
2009). Using these results, we are able to quantify the extent 
to which hydrodynamic bubbles obey simple analytic rela- 
tions and the effect of subgrid turbulence. In particular, the 
inferred sizes of bubbles as a function of the cluster-centric 
distance in our simulations allow us to directly address the 
extent to which such observations can be used to constrain 
the role of magnetic fields in bubble evolution and measure 
the overall energy input from AGN in the centres of cool- 
core clusters. 

The structure of this works is as follows. In §2 we 
describe our simulations, spectral modeling, and analysis 
methods. In §3 we present our results and compare them 
with analytical estimates and observed trends. Conclusions 
are given in §4. 
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2 METHOD 
2.1 Code 

All simulatio ns were performed with FLASH version 3.0 
(jFrvxell et al. 2000), a multidimensional adaptive mesh re- 
finement hydrodynamics code, which solves the Riemann 
problem on a Cartesian grid using a directionally-split 
Piecewise-Parabolic Method (PPM) solver. While the di- 
rect simulation of turbulence is extremely challenging, com- 
putationally exp and dependent on resolution (e.g., 
iGlimm et al.ll200lh . its behavior can be approximated to a 
good degree o f accuracy by adopt i ng a s ub^rid approach. 

Recently, IPimonte Tiptonl (|2006h . described a sub- 
grid model that is especially suited to capturing the 
buoyancy- driven turbulent evolution of AGN bubbles. The 
model captures the self-similar growth of the RT and 
Richtmyer-Meshkov (RM) instabilities by augmenting the 
mean hydrodynamics equations with evolution equations for 
the turbulent kinetic energy per unit mass and the scale 
length of the dominant eddies. The equations are based 
on buoyancy- drag models for RT and RM flows, but con- 
structed with local parameters so that they can be applied to 
multi-dimensional flows with multiple materials. The model 
is self-similar, conserves energy, preserves Galilean invari- 
ance, and works in the presence of shocks. Although it con- 
tains several unknown coefficients, these are determined by 
comparisons with analytic solutions, numerical simulations, 
and experiments. In Scannapieco & Briiggen (2008), we de- 
scribe our implementation of this model into the FLASH 
code, present several comparisons with analytic results, and 
discuss its usefulness for properly capturing the physics of 
AGN-driven cavities in the absence of magnetic fields. The 
reader is referred to this work, for these and other related 
discussions. 

In the present study, all our simulations are performed 
in a cubic three-dimensional region 680 kpc on a side, with 
all reflecting boundaries. For our grid, we chose a block size 
of 8^ zones and an unrefined root grid with 8^ blocks, for 
a native resolution of 10.6 kpc. The refinement criteria are 
the standard density and pressure criteria, and we allow for 
4 levels of refinement beyond the base grid, corresponding 
to a minimum cell size of 0.66 kpc, and an effective grid of 
1024^ zones. 



where r is measured in kpc. Furthermore, the hydrogen num- 
ber density was assumed to be related to the electron num- 
ber density as nu = ne/1.2 according to standard cosmic 
abundances. The static, spherically-symmetric gravitational 
potential was set such that the cluster was in hydrostatic 
equilibrium. 



2.3 Bubble generation 

X-ray cavities in the ICM are thought to be inflated by a 
pair of ambipolar jets from an AGN in the central galaxy 
that inject energy into small regions at their terminal points, 
which expand until they reach pressur e equ ilibrium with the 
surrounding ICM (jBlandford Reeslll974l l. The result is a 
pair of underdense, hot bubbles on opposite sides of the 
cluster centre. 

In order to produce bubbles, we started the simulation 
by injecting a total energy of ^inj into two small spheres 
of radius rinj = 4.5 kpc at distances of 13 kpc from the 
cluster centre. The gas inside these spheres was heated and 
expanded similar to a Sedov explosion to form a pair of 
bubbles in a few Myrs, a time much shorter than the rise 
time of the generated bubbles. The parameters rinj, and 
^inj were chosen such that these regions reached a radius 
of rbbi = 12 kpc and a density contrast of approximately 
Pb/Pamb = 0.05 as compared to the surrounding ICM. How- 
ever, the dependence of the bubble dynamics on the density 
contrast, yOb/pamb, is weak provided that pb/pamb ^ 1- For 
a 7 = 5/3 gas, the energy input during this expansion can 
be simply calculated as 



^expand 



J ^init 



dVp- 



3p evac Vevac 



2/3 



(3) 

where Pevac and 14vac are the pressure and volume at the 
evacuated stage at the end of the expansion and Mnit is the 
initial volume. Finally, the energy released as the bubble 
moves outwards to the edge of the cluster is 



%ou — j 

•^'^init 



dVp- 



3p evac V^vac 



(4) 



The sum of these two contributions for both bubbles 
ergs. 



2.2 Cluster Profile 

For our o verall cluster p r ofile, we adopted the model de- 
scribed in Roedi ger et al ] (|2007l l. which was constructed to 
reproduce the properties of the brightest X-ray cluster A426 
(Perseus), which has been studied extensively with Chan- 
dra and XMM-Newton. In this case, the electron density 
Ue and temperature Te profiles are based on the deprojected 
XMM-Newton data (jChurazov et al.1 l2003l l which are also 
in broad agreement with t he Chandra data (| Schmidt et al.l 
I2OO2I : ISanders et al.ll20o3 ). Namely: 

4.6 X 10~^ 4.8 X 10~^ _3 

and 



2.4 Synthetic Chandra observations and spectral 
modeling 

In order to allow a direct comparison of simulation output 
with X-ray data, we made use of a newly-developed pipeline 
for post-processing of gridded simulation output. T he X - 
ray- imaging pipeline XIM (see also iHeinz Briiggenl l2009l l 
is a publically- available set of scripts that automate the cre- 
ation of simulated X-ray data for a range of satellites. It 
takes as input the density, temperature, and velocity, as well 
as a large number of parameters and provides simulated X- 
ray data in the form of spectral-imaging data cubes. Simple 
batch-processing and trivial parallelisation allow the manip- 
ulation of these large data cubes, provided that the final 
output data fit entirely in memory. 

XIM is focused on visualizing X-ray data from thermal 
plasmas. The scripts allow the choice of a user-supplied 
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spectral model as we ll as the default thermal APEC plasma 
emission model (Sm ith et al.|[20Ql a). which self-consistently 
calculates the equilibrium ionization balance for a ther- 
mal plasma. For a fixed set of abundances and a given 
temperature, APEC then computes interpolated high reso- 
lution model X-ray spectra. Atomic data are taken from the 
ATOMDB using APED (Smith et al. 2001b) and combined with 
bremsstrahlung continuum for all species. 

Querying APED for every simulated cell is computa- 
tionally too expensive. For computational convenience and 
speed, XIM creates a table of model spectra that span a user- 
supplied energy vector and a logarithmic range of temper- 
atures with a user-supplied minimum and maximum. The 
spectral contribution of each cell is then logarithmically in- 
terpolated from this table. In this paper, we have fixed the 
relative abundance of all elements heavier than helium at 
the solar value, but in principle, the metallicity Z of the 
gas can be specified on a cell-by-cell level, for fixed relative 
abundances of heavy elements. 

2.5 Spectral projection 

Having computed the spectral contribution in each cell, XIM 
then calculated a raw spectral data cube by projecting the 
data along one of the three Cartesian coordinate axes. Spec- 
tra were emission-measure weighted and Doppler shifted 
with the user-provided radial velocity, neglecting relativistic 
effects. A user-supplied tracer grid was used to weigh the 
data by the thermal plasma content of the cell. The spectral 
grid was oversampled by a factor of three with respect to the 
final output energy grid to allow accurate representation of 
Doppler shifts in the output spectra. 

The data was further redshifted according to the user- 
specified cosmological redshift and the coordinate axes are 
scaled to the proper angular size, given the redshift and 
cosmological parameters, by default using concordance pa- 
rameters Q = 1, Qa = 0.7, Ho = 70kms~^Mpc~^ The 
flux in a cube (x,y, wavelength) was then scaled to the cos- 
mologically correct flux at the given distance. The pro- 
jected data were processed for foreground photo-electric 
absorption using the Wiscons in Absorption Model (WABS 
iMorrison McCammonlll983l l. 

2.6 Virtual observation 

Next, the raw spectral-imaging cube (x,y, wavelength) was 
re-gridded in the two spatial directions onto the detector 
plate scale of the user-specified instrument. XIM incorpo- 
rates telescope parameters for Chandra, Constellation- 
X, XEUS, and XMM-Newton and wih incorporate a 
telescope model for IXO once response matrices become 
available, in our case the Advanced Imaging Spectrometer 
(ACIS), on board the Chandra X-ray observatory. The re- 
gridded data cube was then convolved with the appropriate 
spectral response and ancillary response matrices, and the 
convolution output was regridded onto the user-specified en- 
ergy grid. 

The output was convolved with a model point-spread 
function for the selected telescope. The current version 
of XIM is limited to a Gaussian point spread function 
with energy-independent kernel width. It takes into ac- 
count quantum efficiency and telescope effective area, but 



neglects detector non-uniformity, vignetting, and point- 
spread function variance. The latter effects can in princi- 
ple be modeled through post-p rocessing with th ird-party 
telescope simulators (e.g., MARX) (jWise et al.lll997l ). Finally, 
Poisson-distributed photon counts were calculated for a 
user-specified exposure time. 



3 RESULTS AND DISCUSSION 
3.1 General Properties 

In Fig. [T] - [2] we show synthetic observations for the ACIS- 
I instrument onboard Chandra that covers the band from 
0.7-10 keV. ACIS is a powerful tool for conducting imaging, 
spectroscopic and temporal studies of X-ray sources. The in- 
strument consists of ten Charge Coupled Devices especially 
designed for efficient X-ray detection and spectroscopy. It is 
the standard instrument to study X-ray cavities. 

We set the cluster to the distance of Perseus (z = 0.017) 
and the observation time to 100 ks, and first considered a 
case in which the initial bubble radii are rbbi = 12 kpc, offset 
from the cluster centre by 13.5 kpc. Figs. [T] and [2] show two 
different angles with respect to the axis of the AGN. In Fig.[T] 
the jet axis, meaning the line connecting the bubble centres, 
lies in the plane of the sky. In Fig. [2l the jet axis is inclined 
by 45 degrees with respect to the line of sight. 

Synthetic observations from a simulation without sub- 
grid turbulence corresponding to Fig.[T]and Fig.[2]are shown 
in Fig. [3] and Fig.|4l respectively. As apparent from the syn- 
thetic observations, the differences produced by the subgrid 
turbulence are significant, especially at times > 150 Myrs. 
Without the subgrid model, a bubble does not form a single 
coherent structure but rather looks patchy, eventually com- 
ing apart into resolution-dependent subclumps as described 
in Scannapieco & Briiggen (2008). The cavities in the run 
without subgrid turbulence also show a weaker X-ray sur- 
face brightness contrast, even though the ambient material 
gets mixed into the bubble by subgrid turbulence. Current 
observations of X-ray cavities support the picture from the 
subgrid simulations, and further observations of bubbles at 
larger distances from cluster centres will help to probe the 
stability of the bubbles and check the predictions from these 
kinds of simulations. 

The turbulence also affects the inferred bubble sizes, 
which have been determined by subtracting a smooth, 
radially-symmetric background model from the X-ray maps, 
as is often done with real observations and as shown in the 
last three panels in Fig. [T] - |4l The corresponding radii are 
shown in Fig. [5]- (6] In most observed cavities the diameter of 
all bubbles is larger than the distance of their centres to the 
AGN, indicating the cavities are relatively young. Hence, in 
our simulations, we only follow their evolution for a time 
up to 205 Myrs, and most of the bubbles that have been 
observed should be compared to the first two panels shown 
in Fig. [T] -|2] At this stage of their evolution, all our bub- 
bles are found to have bright rims. As the bubbles move 
outwards, they quickly become difficult to detect as the X- 
ray emissivity is prop ortional to the square of the density 
(|EnBfin H einz 2002). 

The decrement in surface brightness of a cavity rel- 
ative to the surrounding ICM is a strong function of its 
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Figure 1. Synthetic Chandra observations of a simulated galaxy cluster viewed along an axis perpendicular to the axis of the AGN, 
at times of 50, 100, 155 and 205 Myrs after the launch of the bubble. The first panel shows the synthetic observation. The units on the 
X and 2/- axes are pixels. The greyscale shows the number of counts. The other panels show the X-ray fiux with the image of the initial 
cluster subtracted, which makes it easier to discern the bubbles. 



size and distance from the cluster centre. As described in 
iMcNamara NulsenI (|2007l ). a small bubble of radius R on 
the plane of the sky at a distance r from the cluster centre 
produces a count deficit which goes as R^{1 + T^/rc)~^^, 
where Vc is the core radius of the cluster and (3 a constant. 
The counts from the ICM scale as R^{1 + /r^)-^^^^^^ , 
and the noise in this scales as its square root. Hence, the 
signal-to-noise ratio scales as i?^(l + r^/rc)~^^^^~^^^, which 
means that cavities are easiest to detect when they are large 
and located close to the cluster centre. 



3.2 Evolution of bubble size with radius 

Diehl et al. (2008) provide simple analytic scaling relations 
for bubble size as a function of cluster-centric distance, 
which they used to compare against their compiled observa- 
tional data set. The first is derived from a simple /3-model 
for the gas pressure of the ICM, given by 

_ M 

Pgas(r) =po[l + (r/re)'] \ (5) 

where the central pressure, po, is related to the central den- 
sity, po, and isothermal temperature, T, and average molec- 
ular weight, m, via the equation of state of an ideal gas 
with an adiabatic index F = 5/3. Assuming that a bub- 
ble of initial radius Rb,o expands adiabatically in pressure 
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Figure 2. Synthetic Chandra observations of a simulated galaxy cluster viewed along an axis that is inclined by 45 deg to jet axis. The 
panels show the count rates at times of 50, 100, 155 and 205 Myrs after the launch of the bubbles. 



equilibrium with the surrounding ICM, its radius Rb should 
satisfy 

Pgas(r)i?b^ = Pgas(ro)i?b!o • (6) 

Hence, the bubble radius should evolve with distance from 
the cluster centre, r, as 

Rh{r) = i?,,o[pgas(r)/pgas(ro)]"* (7) 
= Rb,o [l + {r/rcf]^ . 

Thus, at large radii, the cavities should grow asymptotically 
as Rh oc (r/rc)~^^^. For typical values (3 — 0.5 for a 
cluster and F = 5/3 this gives Rh oc (r/rc)^^"^°. 

Conversely, for the bubbles produced by current- 
dominated jets, Diehl et al. (2008) assume a constant cur- 



rent along the interior and exterior of the jet axis, such that 
by Ampere's law B oc IR~^, where / is the current and 
B is the toroidal magnetic field strength. This means that 
p oc oc R~^, which combined with Eq. (|5} means that 
the bubble size evolves as 

Rh{r) - i?6,o[Pgas(r)/pgas(ro)]"* (8) 
= Rb,o [1 + {r/rcf] ' . 

As discussed above, such analytic estimates are not able 
to capture effects such as mass-entrainment and distortion of 
the bubbles by drag forces, which are naturally included in 
the simulation. In fact, our simulations show clearly that the 
cavities evolve aspherically as they rise through the cluster, 
expanding quickly in the perpendicular direction, but ex- 
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Figure 3. Same as Fig. ^ only without a subgrid model for turbulence. 



panding slowly, or even becoming compressed, in the radial 
direction. 

There the radial and perpendicular bubble axes evolve 
very differently as shown in Figs. [5] and [6] The radial extent 
of the bubble changes slowly as a function of distance and 
time, while the perpendicular extent of the bubble changes 
rapidly. In the case where the AGN axis lies in the plane 
of the sky, the semimajor axis of the surface brightness de- 
pression is oriented perpendicular to the gravitational accel- 
eration, and the bubble appears to flatten as it moves out- 
wards. Conversely, in the case where the AGN axis is tilted 
by 45 degrees with respect to the line of sight, the projected 
semimajor axis lies in the projected radial direction, and 
the bubble appears prolate. In both cases, the major axes 
grow much more rapidly than the simple analytic prescrip- 
tion given in Eq. JT]). On the other hand, after about 150 



Myrs, the bubble physically flattens. This becomes visible in 
the projection shown in Fig. [2l Note also, that in this pro- 
jection, the bubble never appears to detach from the centre 
of the image for the time that we have simulated. 

One can attempt to constrain bubble properties by con- 
fronting the theoretical evolution of the bubble size with the 
data. However, we do not try this here and it is difficult 
with the current data. Instead, we merely point out some 
of the complexities associated with such analyses, even in 
the non-magnetic case. In Figs. [5] and [6] we have overplotted 
our measurement s from syntheti c obse rvations with a sam- 
ple of bubbles bv iRaffertv et alJ (|2006l l. All the 64 cavities 
in 32 systems in that sample were observed with Chan- 
dra. The sample ranges in redshift from z =0.0035 to 0.545 
and varies in its composition from groups to rich clusters. 
For each cavity, a size and position were measured, assum- 
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Figure 4. Same as Fig. [2] only without a subgrid model for turbulence. 



ing that the cavity extends to the inner edge of any bright 
surrounding emission. The projected semimajor axis, a, and 
semiminor axis, b, of the cavities were measured by eye from 
the exposure-corrected, unsmoothed images. 

Diehl et al. (2008) infer that the non-magnetic models 
are all too shallow to represent the data accurately, while 
the magnetic models adequately reproduce the steep slope 
of the correlation. However, when we examine the tangen- 
tial and radial sizes separately, it is difficult to discard a 
non-magnetic model. Instead, our simulations show that the 
evolution of the semimajor and semiminor radii is not in 
obvious contradiction to the data. In fact, the slopes that 
we find in Fig. [5] are entirely consistent with the measured 
bubble sizes. In Fig. [6] we have added analytic estimates ac- 
cording to eqs. dll) & dS}, but using the exact pressure profile 
for the cluster in our simulations (p oc nT, with density and 



temperature computed from eqs. [T] & [2|). The size-distance 
relationship of the bubbles is dramatically different from the 
simple analytical estimate of Eq. JT]). Especially the slope of 
the evolution of the semimajor axis is much steeper than the 
analytical estimate for the radius as shown by the dashed 
line. 

In our hydrodynamic model of the cavity evolution, the 
effects of projection are not trivial, and the data points are 
not just systematically shifted toward smaller radii. The 
slopes of the size-distance relationship of the bubbles be- 
come significantly shallower as the AGN axis tilts toward 
the line of sight. The appearance of the bubbles show stark 
differences, especially when the bubbles are further than two 
bubble radii from the centre of the cluster. Fig. [5] suggest 
that most observed bubbles come from systems where the 
AGN axis lies between 90 and 45 degrees from the line of 
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Figure 5. Plot of cavity projected radii versus their projected distance from the cluster centres. The connected squares correspond to the 
measured radii from synthetic observations where the axis of the AGN lies in the plane of the sky. The diamonds correspond to the case 
where the AGN axis is inclined by 45 degrees with respect to the line of sight. The unconnected filled circles give the sizes of observed 
cavities compiled in Rafferty et al. (2006). The bottom row shows the results for a run without a subgrid model. Left: semimajor radii, 
Right: semiminor radii. 



sight. In fact, this is expected geometrically for a random 
distribution of bubble axes. 

In the case where the AGN axis lies at 45 degrees from 
the line of sight, the bubbles look more like horseshoes. 
A confirmation of this will require looking for bubbles at 
greater distances from the centres in clusters where it is be- 
lieved that the jet axis is inclined not too far from the line 
of sight. This could help to verify whether the bubbles are 
mainly hydrodynamic bubbles. One may speculate that the 
giant cavities found in Abel l 2204 resemble thes e late-stage 
horse-shoe shaped bubbles (| Sanders et aljlioosh . 

The spread in bubble radii at a given distance from 
the centre in Fig. [5] - [6] is due to a spread in the initial 



bubble size (or energy) which is determined by properties of 
the AGN and the host cluster. Diehl et al. (2008) find that 
the mass of the supermassive black hole and the core radius 
account for most of the spread. To study the impact of these 
differences further, we carried out two additional runs, both 
of which included subgrid turbulence. In one run, we reduced 
the initial radius of the bubbles to rbbi = 6 kpc, and offset 
them from the centre by 8.5 kpc, and in a second run we 
increased the initial radius of the bubbles to rbbi = 18 kpc, 
offset from the centre by 19.5 kpc. The results of these runs 
are shown in Fig. [T] 

In general, the evolution of the minor and major radii 
are similar across runs. If the AGN axis lies in the plane 
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Figure 6. Same as Fig. [5] only the bubble positions are scaled by the core radius of the cluster rc- The dashed line shows the prediction 
by Eq. (0 for a fixed value of i?b(0) and F = 5/3, and the dot-dashed line shows the prediction by Eq. (|8]). The unconnected filled circles 
give the sizes of observed cavities compiled in Rafferty et al. (2006). Left: semimajor radii, Right: semiminor radii. The bottom row 
shows the results for a run without a subgrid model. 



of the sky, the major axis, which is perpendicular to the 
radial direction, evolves much more rapidly than expected 
from analytic estimates for all choices of initial bubble size. 
Likewise, in this projection, regardless of initial bubble size, 
the minor axis grows only very slowly as the cavity moves 
outward. In all cases, the differences between these runs are 
well within the observational scatter. In the projection where 
the AGN axis is inclined by 45 degrees with respect to the 
line of sight, the major axis, which is now parallel to the ra- 
dial direction, grows at a moderate rate in all three runs. In 
all runs this evolution is much more gradual than measured 
when the bubbles are viewed in the plane of the sky. In all 
runs, the minor axis of the bubbles remain roughly constant 
in this projection, even occasionally shrinking slightly. Thus, 



initial cavity size and radial offset seem the be relatively mi- 
nor issues as compared with the anisotropic evolution of the 
cavities and the angle from which this evolution is viewed. 

Finally, Fig. [8] shows the pV energy as function of pro- 
jected distance from the cluster centre. As discussed above, 
the pV energy inferred from observations appears to increase 
with increasing distance, even though this trend is less pro- 
nounced at small energies. This may in part be due to the 
fact that bubbles tend to be overpressured close to their ori- 
gin, and expand to reach pressure equilibrium as they rise 
though the ICM. However, excess pressure in the observed 
cavities is difficult to measure, except indirectly through the 
presence of shocks and sound waves. 

A second possibility for the radial increase in inferred 



© 2007 RAS, MNRAS 000,[THT3l 



Evolution of X-ray cavities 11 





Figure 7. Semimajor (left) and semiminor (right) projected radii versus bubble positions scaled by the core radius of the cluster Tc- 
The lines correspond to runs with three different sizes for the initial bubbles. Squares: small bubbles (6 kpc), triangles: medium-sized 
bubbles (12 kpc), crosses: large bubbles (18 kpc). The unconnected filled circles give the sizes of observed cavities compiled in Rafferty 
et al. (2006). The top row corresponds to the case where the axis of the AGN lies in the plane of the sky; the bottom row displays results 
for the case where the AGN axis is inclined by 45 degrees with respect to the line of sight. 



pV is the entrainment of ambient material. In fact, this is 
the main reason for the increase in our simulations, as in our 
model the bubbles quickly reach pressure equilibrium with 
the surrounding medium, well before moving noticeably out 
from the cluster centre. Over the 200 Myrs that we simulate, 
we find that the inferred pV energy for the same bubble 
grows by a factor of 2 in the case where the AGN axis lies 
in the plane of the sky. In the case where the AGN axis is 
inclined by 45 degrees with respect to the line of sight, the 
inferred pV energy grows hardly at all. In the latter case, the 
greater apparent pressures (because the bubbles are further 
in) are offset by the bubbles appearing to be smaller. 

The measured energies are more reliable when the bub- 



bles are further from the centre of the cluster, which unfor- 
tunately, is where they are most difficult to observe. Overall, 
the error introduced into estimates of the pV energy by these 
effects is of the same magnitude as the error related to the 
unknown equation of state of the plasma inside the bubbles. 



4 CONCLUSIONS 

The nature of AGN-driven X-ray cavities remains one of the 
major outstanding questions in understanding the physics of 
cool- core galaxy clusters. While the observed morphologies 
and sizes of the cavities provide us with useful clues as to 
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the processes at work in these regions, interpretation of the 
observations is far from straightforward. 

In this study we have focused on a model in which X-ray 
cavities evolve purely hydrodynamically, and reconstructed 
their detailed evolution using of two major tools: AMR sim- 
ulations that include subgrid turbulence and synthetic X-ray 
software that produces realistic observations from these sim- 
ulations. Together these tools allow us to capture such im- 
portant effects as mass-entrainment, distortion of the bub- 
bles by drag forces, and observational effects. These effects 
lead to an evolution that is drastically different than ex- 
pected from simple analytic estimates. 

In particular, we find that while the radial extent of the 
cavities changes slowly as a function of distance and time, 
they expand rapidly in the perpendicular direction. The re- 
sult is a complex evolution that is highly dependent on view- 
ing angle and difficult to compare conclusively with observa- 
tions. Although analytic estimates of non-magnetic models 
evolve too slowly to match the observations, our simulations 
show that the evolution of the semimajor and semiminor 
radii is not in obvious contradiction to the data. In fact, 
our simulations naturally reproduce the overall trend for in- 
ferred pV energies of observed X-ray cavities to increase as 
a function of distance from the cluster centre, an effect that 
is largely due to mixing of entrained material into the rising 
cavities. The size evolution we find is in general good agree- 
ment with the measured bubble sizes, although there is a 
large spread in observational data, and a strong dependence 
on the direction from which the cavities are viewed. Indeed, 
the fiattening and projection of X-ray cavities have much 
stronger effects on our results than initial bubble size and 
radius. 

Finally, some words of caution. Our simulations start 
from a very simplified setup. In reality, the bubbles are in- 
fiated by jets continuously over some time. Also the host 
system will be dynamic, leading to a more complex evolu- 
tion than the simulations presented here can capture. One 
example are the bubbles in Hydra that do not look like the 
bubbles we have shown here. In the Hydra cluster, a rapid 
succession of outbursts from the AGN have continued to 
fuel the bubbles to produce a partially interconnected se- 
ries of bubbles that still have some momentum imparted by 
the original jet. More realistic simulations of specific sys- 
tems that are designed to reproduce all observed features 
are needed in order to make better inferences on the nature 
of X-ray cavities. This, in combination with low- frequency 
radio observations, will help to understand these elusive and 
crucial features of cosmic structure. 
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Figure 8. pV energy as function of projected distance from the 
cluster centre. Again squares correspond to synthetic observations 
where the axis of the AGN lies in the plane of the sky; diamonds 
correspond to the case where the AGN axis is inclined by 45 
degrees with respect to the line of sight. The unconnected filled 
circles give the sizes of observed cavities compiled in Rafferty et 
al. (2006). 
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